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<N ■ ABSTRACT 
> 

We investigate possibilities of solar coronal heating by acoustic waves gener- 
ated not at the photosphere but in the corona, aiming at heating in the mid- to 
low-latitude corona where the low-speed wind is expected to come from. Acous- 
tic waves of period r ~ 100s are triggered by chromospheric reconnection, one 
model of small scale magnetic reconnection events recently proposed by Sturrock. 
These waves having a finite amplitude eventually form shocks to shape sawtooth 
waves (N- waves), and directly heat the surrounding corona by dissipation of their 



wave energy. Outward propagation of the N-waves is treated based on the weak 
shock theory, so that the heating rate can be evaluated consistently with phys- 



ical properties of the background coronal plasma without setting a dissipation 
length in an ad hoc manner. We construct coronal structures from the upper 
chromosphere to the outside of 1AU for various inputs of the acoustic waves hav- 
ing a range of energy flux of F w = (1 — 20) x 10 5 erg cm -2 s _1 and a period of 
t = 60 — 300s. The heating by the N-wave dissipation effectively works in the 
inner corona and we find that the waves of F w ^ > 2 x 10 5 erg cm _2 s~ 1 and r > 60s 
could maintain peak coronal temperature, T max > 10 6 K. The model could also 
reproduce the density profile observed in the streamer region. However, due to 
its short dissipation length, the location of T max is closer to the surface than the 
observation, and the resultant flow velocity of the solar wind is lower than the 
observed profile of the low-speed wind. The cooperations with other heating and 
acceleration sources with the larger dissipation length are inevitable to reproduce 
the real solar corona. 



Subject headings: Sun: corona — solar wind — waves 
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1. Introduction 

The heating of the solar is been still poorly understood, and it is one of the most 
challenging but interesting questions to be solved in astrophysics. The origin of energy that 
heats the corona is generally believed to lurk in the turbulent convective motions beneath the 
photosphere. A certain fraction of the kinetic energy of those turbulent motions is carried 
up to the corona in the shape of non-thermal energy such as magnetic or wave energy, 
and thermalization of such energy in the corona results in the heating of the surrounding 
plasma. Granule motions of the surface convection simply expect wave generation at the 
photospheric level. Possibilities of coronal heating by those waves have been investigated by 
many researchers (Osterbrock 1961; Ulmschneider 1971; McWhirter, Thonemann, & Wilson 
1975 ). To date, among various modes of waves, although an Alfven wave has widely taken 
up as a convincing candidate in coronal heating and acceleration of the high-speed wind in 
the polar coronal holes through the ion-cyclotron damping mechanism (e.g. Cranmer, Field, 
& Kohl 1999; Hollweg 1999), acoustic waves have not been regarded as a major heating 
source of the corona because of their dissipative character. Acoustic waves with a finite 
amplitude inevitably steepen their wave fronts to form shocks when they travel. They are 
strongly damped through upward propagation in the chromosphere, consequently the only 
0.01% of the initial wave energy could reach the corona (Stein & Schwartz 1972; hereafter 
SS). A role of acoustic waves in solar coronal heating has been buried in oblivion for a long 
time. 

However, recent observations of dynamical structures of the solar corona (see Aschwan- 
den, Poland, & Rabin 2001 for recent review) highlight a role of acoustic waves generated 
not at the photosphere but in the corona. A large number of flares and flare-like events (e.g. 
Tsuneta et al.1992; Tsuneta 1996) have been observed by various telescopes until today (e.g. 
Nishio et al.1997; see Bastian, Benz, & Gary 1998 for review) in a wide range of energy from 
~ 10 24 erg to ~ 10 32 erg, and they are statistically argued in terms of supply of bulk energy to 
heat the global corona. Recent measurements of EUV frequency distribution, ^(E) oc E~ a , 
of flare-like events gives a steeper power-law index (a = 2.3 — 2.6 by Krucker & Benz 1998; 
a = 2.0 — 2.6 by Parnell & Jupp 2000) on the lower energy side (E = 10 24 — 10 25 erg). This 
fact makes us infer that small-scale flare-like events, such as nano-flares, might be sufficient 
(Hudson 1991) for the required energy budget (> 10 5,5 erg cm _2 s _1 ; Withbroe & Noyes 1977). 
Several models taking into account such small-scale events have been introduced (Sturrock 
1999; hereafter S99; Sturrock, Roald, & Wolfson 2000; Roald, Sturrock, & Wolfson 2000; 
hereafter RSW; Tarbell, Ryutova, & Covington 1999; Sakai et al. 2000; Furusawa & Sakai 
2000). Among these models, S99 proposed a model of chromospheric magnetic reconnection. 
In this model, a flux tube, newly formed as a result of reconnection events, oscillates verti- 
cally to excite acoustic waves in the corona. RSW further showed based on a simple kinetic 
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model that this process possibly liberates the sufficient energy in the corona with the correct 
relevant parameters. A different mechanism also predicts generation of acoustic waves in the 
corona. For example, spicules in the open magnetic field region effectively transports energy 
of random motions at the photosphere into the corona and produce longitudinal (acoustic) 
waves there (Hollweg 1992; Kudoh & Shibata 1999). Thus, acoustic waves are supposed to 
be generated constantly in the corona by various dynamical processes. Such acoustic waves 
could heat the surrounding corona directly, unlike acoustic waves created at the photospheric 
level. They are expected to heat the inner corona effectively by the shock dissipation. There- 
fore, density at the coronal base might be large enough to creat dense wind flow observed 
as the streamer in the mid- to low-latitude region, which is believed to be connected to the 
low-speed solar wind (Habbal et al. 1997). Acoustic waves might become one of relevant 
processes working in the region generating the low-speed wind of which mechanisms have 
been poorly elucidated. 

In the above models introducing production of acoustic waves in the corona, the authors 
concentrated on the amount of kinetic (wave) energy released in the corona. However, heating 
of coronal gas is accomplished by thermalization of wave energy through dissipation of such 
acoustic waves. Therefore, an appropriate treatment of propagation and dissipation of these 
waves is indispensable in understanding a problem of coronal heating. In this paper, we 
employ a formulation of a weak-shock theory constructed by SS, originally developed to 
study acoustic waves generated by convective motions at the photosphere, for propagation 
and dissipation of the waves produced in the corona. Then, we can explicitly determine the 
rate of heating by dissipation of the waves on given wave energy flux (or amplitude) and 
wave period without tuning other free parameters; we do not have to arbitrarily take an 
exponential type of heating of mechanical energy input, F m oc exp( ~^~ R0 ^ ), for an assumed 
constant dissipation length, / m , which is poorly supported by fundamental physical processes, 
although this conventional shape of the heating law was adopted, for lack of an alternative 
in most previous models (Kopp & Orrall 1976; Withbroe 1988; Sandbaek & Leer 1994), to 
study the global coronal structure. 

The paper is organized as follows: In §2, we present our model. We briefly summarize 
the generation of acoustic waves in the corona (§2.1), and the weak-shock theory formulated 
by SS (§2.2). Basic equations describing the global coronal structure from 1R & to > 1AU are 
shown in §2.3, and we present our method to construct a unique coronal solution on given 
wave parameters. In §3 we show our model results. First, we discuss variation of coronal 
wind structures with respect to input parameters from a theoretical point of view (§3.1). 
Second, we examine several characteristic properties of the resultant corona (§3.2). Finally, 
we test whether our model is able to reproduce recent observed results for the low-speed 
wind in the mid- to low-latitude region (§3.3). In §4, we summarize our results and discuss 
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related topics. 

2. Model 

2.1. Generation of Acoustic Waves in the Corona 

As a new relevant process of coronal heating in the quiet-Sun region, S99 proposed a 
model of network-field magnetic reconnection at the chromospheric level. The chromosphere, 
specially at the minimum-temperature location, is quite a favorable site for reconnection to 
occur, since the magnetic resistivity is greatest, and the width of the current sheet, which 
is possibly scaled with the pressure scale height, is smallest there. Such reconnection events 
form a new closed magnetic flux tube, which is initially located far from the equilibrium 
state, and hence, will spring upward quickly (consult S99 for a schematic picture). The 
tube would oscillate about the equilibrium state normal to the solar surface with a period of 
r = 2L/va, where va is Alfven velocity and L is the length of the tube. If we use the quoted 
values in S99 as a typical example of the tube, L ~ 1.5 x 10 9 cm, a mean magnetic field, 
~ 100G, and a mean density, ~ 10 _12 g cm -3 , the oscillation periods would be r ~ 100s. 
It could keep oscillating several tens of times before it is damped (S99). As a result, these 
perpendicular oscillations would excite longitudinal waves traveling in the vertical direction. 
Given that the configurations of the magnetic fields are perpendicular above the flux tube, 
such waves would propagate upwardly as acoustic waves in the rarefied atmosphere. They 
must contribute to the heating of the surrounding corona directly by the dissipation of the 
waves. 

Using a simple kinetic model, RSW have estimated the energy liberated by the above 
mechanism as a function of a mean magnetic field strength at the photospheric level. They 
showed that the mean field of ~ 10G could generate the required heating rate (> 10 5 5 erg 
cm _2 s _1 ) to explain the conditions in the quiet-Sun region. A sizable fraction of the released 
energy is supposed to be transported to the energy of the acoustic waves. In our model, 
we parameterize the input energy flux of the acoustic waves (not the total energy liberated 
by the reconnection) as F w0 (erg cm _2 s _1 ). In this paper, to study the physical processes 
clearly, we focus on the role of such acoustic waves in the coronal heating and construct 
the coronal structure on the input F Wj0 and wave period r, though some of the released 
energy would be transformed into other MHD waves depending on complex configurations 
of magnetic fields (Tarbell, Ryutova, & Covington 1999; Sakai et al. 2000). 

We would like to remark that other mechanisms also predict production of longitudinal 
waves at the coronal height, although we have taken up the chromospheric reconnection 
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model as a typical process in this paper. Generation of spicules in an open magnetic flux 
tube have been investigated by various authors (Hollweg, Jackson, & Galloway 1982; Hollweg 
1992; Kudoh & Shibata 1999). They found that Alfven waves excited by random motions at 
the photosphere (Ulrich 1996) effectively transport their energy to the corona. The nonlinear 
effect of torsional Alfven waves produces longitudinal waves along the vertical flux tube in 
the corona. A sizable fraction of the initial energy of the transverse waves at the photosphere 
is converted to energy of the longitudinal waves at the coronal height, and these waves could 
become acoustic waves propagating upwardly. They could heat the surrounding plasma in 
the very same way as those triggered by the reconnection events above. Thus, acoustic waves 
are expected to be universally generated by various mechanisms in the corona far above the 
photosphere, and therefore, it is quite worth studying their role in coronal heating. 



In this section, we describe our method of treating the outward propagation of acoustic 
waves, after they have been excited in the corona. We first estimate a distance acoustic waves 
travel before forming shocks on a plane-parallel geometry. Then, we derive an equation 
dealing with variation of N-wave amplitude on spherical geometry. For simplicity's sake, 
in the following discussions, we neglect the effects of magnetic field on the propagation of 
the waves. This simplification is valid when the circumstantial magnetic configuration is 
perpendicular above the flux tube generating the waves. 



Any acoustic wave having a finite amplitude inevitably changes its shape, makes the 
wave front steepen and eventually forms the shock front (e.g. Landau & Lifshitz 1959). The 
distance the acoustic waves travel before forming the shocks can be estimated on a given 
wave length, A, and initial amplitude, 5vq. Consider the acoustic waves propagating in the 
upward direction, z, in isothermal atmosphere with density structure of p = poexp(—z/H p ). 
Assuming waves having initially sinusoidal velocity profiles, 5v = 5v sm(2irZ/ \), the wave 
crest overtakes the preceding trough to form a shock front at 



(SS), where z is the position at which the waves are created, 7 = 5/3 is a ratio of specific 
heat, and c s is a sound velocity. 



2.2. Dissipation of Acoustic Waves in Corona 



2.2.1. A Distance Acoustic Waves Travel Before Forming Shocks 




(1) 
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In our calculations, waves of the initial amplitude, Svo/c s = 0.1 ~ 1 are considered 
(§3.1), and therefore, the second term in the logarithm in eq.(l) will be bounded by an 
upper limit: 

4(7 + l)^5t ~ W p = ( ^ rl ° sr ~ °- 25( 100i )( 2xloims-i rl (2) 

where g is acceleration of gravity, and we have used H p ~ H p which is satisfied in the corona 
where temperature varies slowly within a scale of A(~ 10 9 cm). Then, eq.(l) can be expanded 
to a first order as 

A c s 

2 " 20 = 2(7TT) cV (3) 

The factor 2 ( 7 +i) * s an order of unity, indicating that after traveling one wave length, the 
waves begin to dissipate energy. 



2.2.2. Variation of N-wave amplitude 

The initially sinusoidal acoustic waves are expected to be transported as N-waves after 
the formation of the shocks, provided that they are continuously generated from the lower 
corona. These N-waves are propagated outwardly in the rarefied atmosphere and dissipate 
their wave energy to heat the corona. For the purpose of giving a reasonable estimate of 
heating rate as a function of position, variations of amplitude, 5v w , of the transported N- 
waves are treated based on the weak-shock theory, following the formulation presented by 
SS. In the discussions below, propagation of the wave train are considered in a flow tube with 
a cross-section of A, to keep the consistencies with a model for the global corona presented 
in §2.3. A is a function of r, a distance measured from the center of the Sun, and modeled 
in §2.3 with taking into account non-radial expansion of the flow tube. From now on, all the 
physical quantities are expressed as functions of r only, unless explicitly declared. 

An equation for the variation of wave amplitude normalized by ambient sound velocity, 
a w = — , can be found from SS: 



1 da w 1 . 1 dp 1 dE\ 1 d\ . . 4 , 

+ 77- -T^-T-y). (4) 



a w dr 2 p dr E\ dr X dr 

where p is gas pressure, and E\ = |p(<5t> w ) 2 A = ^jpa^X is a wave energy per wave length, 
A. The variation of the wave energy can be estimated from entropy generation by the weak 
shock (SS; §56 in Mihalas & Mihalas 1984) as 

dE\ E\ dA 2( 7 + IK 2( 7 + l)q w 

V ' Ex = ~dV + -A^ = - Ex A " cr ' (5) 
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where we have used a relation of A = c s (l + : ^-a w )r ~ c s r by assuming a w < 1. In general, 
the period of waves traveling in different media remains as a constant, implying = ^-^f- 
Hence, eq.(4) is reduced to 

da w _ a w I dp 2(7 + l)a w 1 dA 1 dc s 
dr 2 p dr c s r A dr c s dr 

This equation determines the variation of the wave amplitude in the solar corona according 
to the physical properties of the background coronal plasma. The first term in the r.h.s. is 
positive in the density decreasing atmosphere, and the second term of the entropy generation 
(heating ) is negative. In the (1 — 2)R Q region, where the dissipation is important to heat the 
corona, these two terms always dominate the other terms, the third term arising from the 
geometrical expansion and the fourth due to temperature variation. Using a w determined 
above as well as background physical quantities, p and c s , wave energy flux, F w (erg cm~ 2 s _1 ), 
is derived as 

F w = -p C 3(a w ) 2 (l + l±-a w ), (7) 

if recalling that the wave crest moves at a speed of c s (l + ^a w ). 

Since the model described above is simple, we instead have several limitations that 
should be taken into account with great care. The first obvious limitation is that the wave 
amplitude should satisfy the assumption of weak-shock, a w < 1. Second, it does not take 
into account the effects of gravity. Third, it is constructed in a static medium, therefore we 
cannot apply it to cases of moving media, such as N-waves in the solar wind where the flow 
velocity, v, exceeds the ambient c s . Of the three limitations, the third one is most easily 
overcome, because in all the cases we calculate in this paper, at least 99% of the initial wave 
energy dissipates within r < 1.3R Q , the region where v <C c s is fulfilled (figs.l & 2). The 
second limitation also seems to have little effect if the wave period is small enough, since SS 
found that the weak-shock theory gives reasonable estimates for waves of r < where 
c<j ac = is the acoustic cut-off frequency in a gravitationally stratified atmosphere, after 
comparing results of the weak-shock theory with those of fully non-linear calculations. As 
c<j ac ~ 80 o( s ) ^xic/cms- 1 ) _1 m ^ ne corona, the weak-shock theory seems to be applicable to 
waves of r < 2000s. To check whether the first limitation is overcome, we need to calculate a 
a w variation in the solar corona at first hand. Our results show that any waves we calculate 
give a w < 0.5 in the entire corona (§3.1). Thus, the weak-shock theory appears to be 
applicable to our model and to give a reasonable estimate of the heating rate as a function 
of r. 
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2.3. Basic Equations 

We here present basic equations to describe one-component coronal wind structure in 
a flow tube with a cross-section of A under a steady state condition. Then, an equation of 
continuity becomes 

pvA = const. (8) 
An equation of momentum conservation is 

dr pdr pc s (l + I±i« w ) V ^ W 

Pressure, p, is related to p and temperature, T, by an equation of state for ideal gas: 

P = P—T, (10) 

where mn is hydrogen mass, k B is Boltzmann constant, and p is mean atomic weight of 
particles in unit of mu- The third term in eq.(9) represents wave momentum deposition, 
neglecting the effects of reflection and refraction (Rosner & Vaiana 1977). Total energy 
equation is obtained as 

V • \p»{\f + -J_^T - ^® ) + F w + F c ] +q R = 0, (11) 
2 7 — 1 m H /i r 

where we adopt the classical form of conductive flux for ionized gas, 

dT 5 dT 

Fc = _ K _ = _ Ko r,- , (12) 

with /to = 1.0 x 10~ 6 in c.g.s unit (Allen 1973), and the radiative cooling term, (/R,(erg 
cm~ 3 s _1 ), is derived from tabulated radiative loss function, A(erg cm 3 s _1 ), (Landini & 
Monsignori-Fossi 1990) for the optically thin plasma as g R = n e n p A, where n e & n p are 
electron and H + density respectively. They are calculated by solving ionization of H and He 
(we ignore heavy elements) under LTE condition. The term, V-F w (< 0), in eq.(ll) indicates 
the heating by the dissipation of the waves, which is explicitly written as 

JO r<r d 
V - Fw -\p^(l + ^aJf(f) + (fa w + ^)(f)^] r>r d ' ^ 

where r d is the position where the waves start to dissipate by shaping N-waves, and the height 
hd = rd — Rq from the photosphere corresponds to the sum of the height of the magnetic 
flux tube (~ 10 9 cm; the generation point of the waves) and the distance the acoustic waves 
travel before forming the shock fronts (~ A ~ 10 9 cm; eq.(3)). We adopt h d = 2 x 10 9 cm as 
a standard value. 
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To take into account non-radial expansion of the flow tube due to configurations of 
magnetic field, cross-sectional area, A, is modeled as 

f P (r-ri)/(T — I — f 
A = r 2 /maxe ±Jl (14) 

where 

/l = 1 - (/max " 1 y 1 -^ 

(Kopp & Orrall 1976; Withbroe 1988). The cross section expands from unity to / max most 
drastically between r = r\ — a and r\ + er. Of three input parameters, / max is the most 
important in determining solar wind structure. In this paper, we consider cases for / max = 5 
and / ma x = 1 (pure radial expansion). As for the rests of the two parameters, we employ 
r\ = 1.25Rq and a = 0.1R Q , the same values adopted in a model for 'quiet corona' in 
Withbroe (1988). 



2.4. Boundary Conditions and Computational Method 

Now we would like to explain the practical aspects of our method of constructing wind 
structure with respect to various input properties of the waves. In order to solve both 
heating of corona (energy transfer) and formation of solar wind (momentum transfer) con- 
sistently, our calculation is performed in a broad region from an inner boundary of the 
upper chromosphere where temperature, = 10 4 K, at r c h = R@ + h^, which is located 
at /i c h = 2 x 10 8 cm (2000 km) above the photosphere (Allen 1973) to an arbitrary outer 
boundary at r out = 3OO-R . 

We set four boundary conditions to construct a unique solution for a given input wave 
energy flux F w as follows: 

F w (r ch ) = F w (r d ) = F w , , (15) 

T(r ch ) = T ch , (16) 

|F c (r ch )|(~0) < |F c , max |, (17) 

V • F c (r out ) = 0, (18) 

where -F c ,max in e Q-(17) is a maximum value of downward conductive flux in the inner corona. 
The first condition denotes that wave energy flux must agree with a given value when the 
waves start to dissipate. The second condition is also straightforward: the temperature 
has to coincide with the fixed value at the inner boundary. The third condition is the 
requirement that the downward thermal conductive flux should become sufficiently small at 
the upper chromosphere (T = 10 4 K), diminishing from its enormous value at the coronal 
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base (T ~ 10 6 K). Practically, we continue calculations iteratively until F c (r c t) / F CjTaax < 1% 
is satisfied. The fourth condition corresponds to an ordinary requirement that no heat is 
conducted inward from infinity (Sandbaek & Leer 1994). Note that thanks to the third 
condition, coronal base density, which is poorly determined from the observations, does not 
have to be used as a boundary condition. As a result, the number of free parameters to be 
set in advance is reduced (Hammer 1982a,b; Withbroe 1988). The density at the coronal 
base or the transition region (TR) is calculated as an output; larger input F w increases 
downward F c in the lower corona, demanding larger density in the coronal base and TR to 
enhance radiative cooling to balance with the increased conductive heating. 

For numerical integration of the momentum equation (eq.(9)) and the energy equation 
(eq.(ll)), we respectively use v and an isothermal sound velocity, a, defined as 

a 2 = °1 = V - = J*-T. (19) 
7 P map 

To carry out the integration, the equations shown in the previous section need to be trans- 
formed into useful forms. First, an expression for velocity gradient can be written from 
eqs.(8),(9), & (10): 



An, _GMq , a?dA _ da? 
_ ~? r ~ ~r A dr dr 

dr v- — 

v 



(20) 



Second, an expression for gradient of a 2 is derived from an integrated form of eq.(ll): 

= ^(^)H V 4 + - Jy - ^ + ^ + f dr<* - £„) , (21) 
dr K a° m^p 2 7 — 1 r pv J ro pv 

where, r is a certain reference point of integration that will be set later. E tot is mathemat- 
ically an integral constant, which must be conserved in the entire region of the calculation. 

Only transonic solutions are allowed for the flow speed, v. The integration of three 
differential equations (20), (21), and (6) is carried out simultaneously from the sonic point, 
r = r s to both outward and inward directions by the fourth-order Runge-Kutta method, also 
deriving density from eq.(8) at each integration. We have used variable sizes for the grid 
of the integration, setting the smaller mesh size in the region where physical values change 
rapidly. For instance, in the TR, we set a mesh as small as 10 3 cm (0.01 km) for one grid. 
To start the integration, we have to set nine variables, v s , (^) s , a 2 , (^§r) s , tt WjS , (^f^) s , Ps, E tot 
and r s (subscript 's' denotes the sonic point). We can determine v s , (^) s and r s for given a 2 
and (^f-)s by a condition that both the numerator and the denominator of eq.(20) are zero 
at r s (Parker 1958). (^)s is also derived from eq.(6) on a given a w>s . Setting the reference 
point, r = r s , for the integration of radiative cooling function, we obtain E tot as 

r l o 7 o GM G) F w ,mup.r 1 ^da 2 . , 

^tat = hv 2 + —-a 2 ^ + — - )*-K a 5 (-r- s , 22 

2 7 — 1 r pv k B pv dr 
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where all the variables are evaluated at the sonic point. Now we have four variables, 
°s> (t~)s) a w,si an d Psj remaining to be regulated by the four boundary conditions of eqs.(15)- 
(18). Concrete procedures for finding a unique solution on a given F w are described below. 

1. One makes an initial guess for [a%, (^r) s , a w,s, As]- 

2. The integration is performed in the outward direction from r s . Leaving a 2 ,a WiS ,&p s 
unchanged, (^§r) s are determined by carrying out the integration iteratively to satisfy 
the outer boundary condition of eq.(18). 

3. The integration is performed iteratively in the inward direction, improving a 2 , a W)S , 
and p s for the fixed (^§r) s until they satisfy the three inner boundary conditions, 
eqs.(15)-(17). Physically, the condition of the wave flux (eq.(15)) regulates a WjS , that 
of temperature (eq.(16)) regulates a 2 , and that of the conductive flux (eq.(17)) regulates 
p s . These relations guide improvement of the respective initial guesses, though they are 
not independently approved. Unless the above inner boundary conditions are satisfied 
simultaneously, one returns to procedure 2, preparing a new set of [a 2 , (^-) S) tt WjS , p s ]. 

4. One can finally find a unique solution by iterating procedures 2 and 3. 

3. Results 

3.1. Theoretical Interpretation of Resultant Coronal Structures 

Figure 1 demonstrates our result of coronal wind structures employing the same input 
-^w,o — 7.8 x 10 5 erg cm -2 s -1 , but different sets of wave periods and non-radial expansion fac- 
tors, (t(s), /max) = (300, 5), (60, 5), (300, 1). Figure 2 compares the wind structures adopting 
three different F w = (3.2, 5.9, 10) x 10 5 erg cm~ 2 s _1 for identical inputs of r = 120s and 
/max — 5. Distribution of temperature, flow velocity, and electron density from the inner 
boundary to 1AU (= 215i? ) are displayed from top to bottom on the left side of the both 
figures. On the right side, we present variation of wave amplitude, a w , a dissipation length, 
l w , of N-waves, which is defined as l w = ^fp-, and heating per unit mass, |^V • F w |, respec- 
tively. In tab.l, we tabulated several resultant properties of the corona and solar wind as 
well as the input parameters. 

To begin with, we would like to emphasize that corona heats up to > 10 6 K in every 
case, because acoustic waves generated in the corona are able to heat the surrounding gas 
directly, unlike acoustic waves produced at the photosphere. However, N-waves are rapidly 
damped, so that the heating occurs only in inner region as seen in the lower right panels of 
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Input 


Out 


put 


\ r vffii 1 j J max) 


Ptr 


T 

± max 


r Tmax 


(n p w)iAu 


(u)lAU 


(7.8, 60, 5) 


0.28 


1.43 


1.04 


3.8 x 10 6 


222 


(7.8, 300, 5) 


0.24 


1.50 


1.08 


7.5 x 10 7 


251 


(7.8, 300, 1) 


0.29 


1.71 


1.12 


3.8 x 10 8 


222 


(3.2, 120, 5) 


0.13 


1.16 


1.05 


2.5 x 10 6 


230 


(5.9, 120, 5) 


0.21 


1.35 


1.05 


1.3 x 10 7 


240 


(10, 120, 5) 


0.32 


1.52 


1.05 


1.3 x 10 7 


230 



Table 1: Input parameters and output wind properties of each model; F Wy o is in 10 5 erg 
cm _2 s _1 and r in second. p tr is pressure in dyn cm -2 at the TR where T = 10 5 K, T max is 
peak coronal temperature in 10 6 K, r Tmax is location of T max in R Q , (n p t>)!Au is proton flux 
in cm _2 s _1 at 1AU, and (f)iAu is the flow velocity in km/s at 1AU. 

figs.l and 2. Then, location of maximum temperature, T max , is quite close to the surface. In 
an outside region of > 1.5i? Q , heat is input only by outward thermal conduction, and the 
flow is accelerated mostly by thermal pressure. As a result, speed of the solar wind at 1AU 
is < 300km/s, which is slightly slower than the actual low-speed wind (300 ~ 450km/s). 

The top right panel of fig. 2 interestingly illustrates that distributions of a w are almost 
identical in spite of very different inputs of F w . Particularly, initial N-wave amplitudes, 
a w (r d ), at r d , are within a range between 0.48 and 0.49. This is because F Wj0 (~ pa^ c s ; 
eq.(7)) mostly owes its variation to change of ambient pressure (see §3.2.1). Moreover, an 
upper right panel of fig. 1 also indicates that initial a w (~ 0.5) is almost independent of 
t and /max- We have found that 0.45 < « w (^d) < 0.52 within our parameter regions of 
F w fl = (1 — 20) x 10 5 erg cm _2 s _1 , r = 60 — 300s, and / max = 1 — 5. This proves that a w < 1 
is fulfilled in the entire region, which justifies the assumption of weak shock. Middle right 
panels of figs.l and 2 show that the dissipation length is a drastically varying function on r. 
This implies that the assumption of a constant dissipation length usually taken in previous 
models for the global corona (Withbroe 1988; Sandbaek & Leer 1994) is very poor for our 
N-wave process. 

In the following discussions, we examine dependences of the wind structures on the re- 
spective input parameters. First, we argue dependences on wave periods. As illustrated in 
fig. 1, N- waves with smaller r dissipate more quickly and the heating occurs in thinner region 
close to the surface. This simply leads to deposition of wave energy in denser region. Since 
radiative loss, (/R(erg cm _3 s _1 ), is in proportion to p 2 for optically thin plasma, a greater frac- 
tion of energy supplied in denser region goes into radiative escape. Consequently, a smaller 
amount of energy remains to heat the corona and accelerate the flow. The case adopting 
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smaller r(=60s) gives lower temperature in the corona, and therefore, a smaller pressure 
scale height and a more rapid decrease of density, as shown in fig. 1. Lower temperature also 
takes the sonic point more distant from the solar surface, and then, mass flux of the solar 
wind becomes much smaller than that expected from the r = 300s case (tab.l). 

Second, we study effects on areal expansion of the flow tube. Comparing results adopting 
the same F w0 = 7.8 x 10 5 erg cm _2 s _1 and r = 300s but different / max = 1 and 5 in fig. 1, 
one can notice significant change of density structure. The model considering the non-radial 
expansion gives more drastic decrease of density as a function of r in spite of similar initial 
density at the inner boundary. Temperature in the inner corona is also lower in that model, 
since more fraction of the input energy is lost adiabatically due to geometrical expansion of 
the flow tube. On the other hand, decrease of temperature is slower and temperature in the 
outer region (> 8R G ) is higher. This is because the lower density reduces both radiative 
cooling and adiabatic loss. The higher temperature also leads to larger acceleration of the 
flow there, giving larger speed of the solar wind in the outer corona. 

Third, we investigate dependencies on input wave energy flux. Models employing larger 
input F w give larger density in the inner corona (fig.2), because radiative cooling should 
be raised to offset the increased heating. More accurately, larger input F Wi0 results in larger 
downward conductive flux from the corona to the chromosphere, which needs larger radia- 
tive cooling to balance with enhanced conductive heating. Larger F w o also leads to higher 
temperature in the inner corona. However, more rapid decrease of temperature occurs be- 
cause of enhanced radiative escape. A model employing F Wi0 = 10 x 10 5 erg cm _2 s _1 yields 
lower temperature in a region of r > 1.5R Q than that adopting F Wi0 = 5.9 x 10 5 erg cm _2 s _1 . 
Therefore, a larger input of F Wj0 does not simply anticipate larger acceleration of the flow, 
and relations of F w — (^pw)iAu an d F w — (i>)iau does not show simple positive correlations 
(tab.l, see also §3.2.2) 

3.2. Some Characteristic Properties of Resultant Corona 

3.2.1. Peak Temperature and Pressure at the TR 

We would like to inspect relations of peak coronal temperature, T max , and pressure, p tT , 
at the TR with respect to the input model parameters. In fig. 3, we present T max as a function 
of input F W) o for different r and /max- The figure shows that T max > 10 6 K is accomplished 
for F w > 2 x 10 5 erg cm _2 s _1 , even if one chooses waves with a short period (r = 60s) and 
a large expansion factor (/ max = 5). T max is a monotonically increasing function of F w on 
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each set of (r, / max ) , and the relations can be fitted by power-law as 

T max oc {F wfi ) k , where k = 0.23 - 0.26. (23) 

Figure 4 displays relations between F Wj0 and p tI . p tI is evaluated at T = 10 5 K, whereas 
p tI is almost a constant through the TR from the upper chromosphere (T ~ 10 4 K; inner 
boundary) to the coronal base (T ~ 5 x 10 5 K) because of its geometrically thin configuration. 
The figure indicates that p tr weakly depends on r and /max, and has a positive correlation 
with F wfi : 

p tT oc (F Wi0 )', where / = 0.75 - 0.79 (24) 

A rise of F Wj o leads directly to an increase in downward conductive flux, which demands 
higher density (or pressure) in the lower corona and the TR to raise radiative cooling to 
balance with the enhanced conductive heating. Alternatively, it could be interpreted that 
injection of larger energy in the corona can heat more deeply to the higher density chromo- 
sphere by the downward thermal conduction. 

A combination of the above two relations of eqs.(23) and (24) roughly give 

T max oc (p tr ) 1/3 , 

on given r and f max , which reminds us of the famous RTV scaling law (Rosner, Tucker, 
& Vaiana 1978), T max ~ 1400(piL h ) 1 / 3 , for closed magnetic loops, where pi(dyn cur 2 ) is 
loop pressure and L h (cm) is loop height. Taking r Tmax instead of L h for our model, we can 
actually derive a relation between T max and (ptr^Tmax) from our results as 

T max ~ 3000(p tr r Tmax )°- 30 , 

showing a form analogous to that of the original RTV law. This is because we consider 
the same energy balance among thermal conduction, radiative cooling, and heating with the 
same boundary condition that the conductive flux should become almost zero at the base 
(eq.(17)), although the configurations are quite different (closed loop for the RTV law and 
open flow tube for ours). The slight discrepancies of the prefactor and power-law index are 
caused by the fact that the RTV law was derived on the assumption of spatially uniform 
heating along the loop, while our heating function is determined by eq.(13), which is not 
uniform at all. 



3.2.2. Mass Flux and Coronal Energy Loss 

In fig. 5 we show anticipated proton flux, (n p i>)iAu, at 1AU as a function of F Wj0 for 
different r and / max , with observational constraints, (n p v )iau = (3-8 ± 1.5) x 10 8 cm~ 2 s~ 1 
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(Shaded), compiled by Withbroe (1988) as the empirical value for 'quiet corona' which is 
supposed to corresponds to the mid- to low-latitude region generating the low-speed wind. 
Larger r waves give a greater (u p v)iav owing to the effective transport of dissipated energy 
to solar wind flow by avoiding radiative escape. Introduction of the non-radial expansion 
of the flow tube reduces (n p v)\Av, because input energy per unit flow tube normalized 
at 1AU decreases on increasing / max , even though one inputs identical wave energy flux 
at the inner corona. Therefore, the larger areal expansion straightforwardly reduces mass 
flux of the solar wind. As to dependences on F w , (n p v)i\u has an upper limit for given 
r and / max - (Even models of r = 300s and / max = 1 are supposed to have the upper 
limit in a F Wj0 > 2 x 10 6 erg cm" 2 s _1 area.) Figure 5 shows that only one case employing 
(ji /max) = (300, 1) can reproduce observed proton flux of the slow wind. Unfortunately 
more realistic cases considering the non-radial expansion of / max = 5 cannot explain the 
observations, which implies that other mechanisms of heating (and acceleration) necessarily 
work cooperatively. 

We would like to study in detail why (n p v )iav has an upper limit on fixed r and / max - 
Figure 6 displays the ratios of the main three types of energy loss: downward thermal con- 
duction from the coronal base, radiative escape in the corona, and total amount of energy 
converted to the flow (mass loss) at r out , respectively normalized by F w for the cases adopt- 
ing r = 120s and / max = 5. The values for downward thermal conduction in fig. 6 are taken 
from maximum conductive flux F c max in lower corona (see fig. 7) and the values for the radia- 
tive escape are derived by the integration of the radiative cooling function from points with 
-F c ,max to r out . The region not labeled between 'radiative escape' and 'flow' denotes energy 
carried out of r out by thermal conduction. According to fig. 6, the main source of coronal en- 
ergy loss is downward thermal conduction within our range of F W) q (strictly speaking, most 
of the conducted energy finally radiates away), whereas radiative escape comes to play a 
significant role for larger F w since density in the corona becomes higher, being subject to 
eq.(24). Although a ratio of energy transfered to the flow is as small as < 3% throughout the 
range, it has a bimodal tendency. In F Wj0 < 6 x 10 5 erg cm~ 2 s~ 1 it increases, which implies 
that the larger F Wi0 is, the more effectively the energy is transfered to the flow. However, 
in F W)0 > 6 x 10 5 erg cm" 2 s _1 , it decreases because of the abrupt dominance of radiative 
cooling. As a result, the predicted (n p v )iau increases rapidly on increasing F w at first and 
eventually decreases as seen in fig. 5. 

To examine these differences in terms of energy transfer, we show variations of energy 
flux of four components, 



wave 




(25) 
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A(r) 

conduction : f c = — ^ F c , (26) 



fl f A^L\( l 2^7 2 GM Q 1 2 7 2 GM Q 
flow : /f = ^^[(^ + " ~ } " + ^1° " ~ ^ ^ 



radiation : / R = pv— — - / dr — , (28) 
A(r ch ) J Tch pv 

per flow tube with a cross section of A =lcm 2 at the inner boundary. The 'flow' term of 
eq.(27) contains three ingredients, kinetic energy of the solar wind, enthalpy, and gravita- 
tional energy. Note that, except the wave energy flux, the zero point of the energy is taken 
here at the inner boundary, to clarify what fraction of the input energy flux is transfered 
to the other components. In a left panel of fig. 7, the status of energy transfer in the inner 
corona is displayed in a linear scale for both X and Y-axis, and in a right panel, that in 
the broader region is shown in log scale. The ratio of the downward thermal conduction, 
f c /F w fi, is smaller for larger F w , though the absolute value of F c is increasing along with 
F w . At the TR, most of the heat flux by the downward conduction finally escapes as radi- 
ation, except for a tiny fraction transfered to the enthalpy to be used to heat the TR. This 
leads to a smaller ratio of the radiative loss, /r/-F w ,o, for larger in the inner region of 
< 1.2i? Q . However, the contribution from the radiation continues in a much more distant 
region in the model employing largest F Wj0 , and /r/-F w ,o finally outdoes the other two cases. 
Consequently, the ratio of energy transfered to the solar wind becomes smaller than the 
model adopting smaller -F w ,o(= 5.9 x 10 5 erg cm' 2 s _1 ) as seen in the right panel. It can 
be concluded that if F w is larger than a certain threshold, an increase of the input F w 
does not lead to effective heating of the corona to accelerate the solar wind but results in 
deposition of the wave energy in the high density region to be wasted as radiative escape. 



3.2.3. Wave Amplitude in the Inner Corona 

Ultraviolet and X-ray emission lines of the corona show non-thermal broadenings (Has- 
sler et al. 1990; Erdelyi et al. 1998), which are inferred to originate from wave motions. We 
investigate whether the obtained wave amplitudes are consistent with these observations. In 
tab. 2, we show our results of the wave rms velocities in the inner corona, which are calculated 

as 



< <5 W >= 




for the N-shaped wave. As the waves propagate upwardly from a location of T = 5 x 10 5 K to 
that of T = 10 6 K, the wave amplitude increases according as the ambient density decreases. 
Our models give the very similar results for different F Wi0 (also r and /max), since change of 
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Input 


<s w > 


\ r vffii 1 j J max) 


T = 5 x 10 5 K 


T = 10 6 K 


(7.8, 60, 5) 


32 


42 


(7.8, 300, 5) 


34 


44 


(7.8, 300, 1) 


32 


41 


(3.2, 120, 5) 


31 


41 


(5.9, 120, 5) 


32 


42 


(10, 120, 5) 


34 


44 



Table 2: < 5 W >(km/s) in the inner corona 



-^ w ,o(~ p < S w > 2 Cg) mostly reflects variation of the density (fig. 2 and §3.1). The tabulated 
results are marginally consistent with non-thermal velocities of 20 ~ 40km/s obtained from 
observations of inner coronal lines in the solar disks and limbs (Erdelyi et al. 1998), if we 
assume the observed non-thermal components totally consist of the waves of the acoustic 
mode. However, this assumption may be too extreme, because other modes of the waves 
actually exist in the real solar corona. 

3.3. Comparison of Coronal Wind Structure with Observed Low-Speed Wind 

In this section, we study a feasibility of the process of acoustic waves in the coronal 
heating by comparing our results with recent observations. Aiming at coronal heating and 
wind acceleration in a region where the low-speed wind is formed, we take observational data 
of the coronal streamer in the mid- to low-latitude region. 

3.3.1. Density Distribution 

In fig. 8, we display our results of electron density with observation of the streamer 
(Parenti et al. 2000; Hayes, Vourlidas, & Howard 2001). Our models adopt an identical 
F Wj0 = 7.8 x 10 5 erg cm'V 1 and three sets of (r(s),/ max ) = (300, 5), (60, 5), (300, 1) (same 
as fig.l). As for the observation, we show results derived from a line ratio of Si IX by 
CDS/SOHO for a region of 1.02 — 1.19i? Q , a ratio of radiative and collisional intensities of 
O VI line by UVCS/SOHO for a region of 1.58 - 1.6-R© (Parenti et al. 2000), and total 
brightness obtained from LASCO/SOHO for a region 1.5i? < r < 6R & (Hayes, Vourlidas, 
& Howard 2001). With respect to the CDS and UVCS data, both results for equatorial and 
mid-latitude streamer regions are displayed. 
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Our results, adopting the same F w $, show almost identical density in very inner part, 
being independent of r and / max (§3.2.1 and fig. 4). The figure exhibits that the adopted value, 
F w = 7.8 x 10 5 erg cm~ 2 s _1 , reproduces the CDS observation of the mid-latitude streamer 
well. To fit to the data of the equatorial streamer, slightly smaller input of F w is favored. 
In our model, r and / max control slope of density distribution, r influences temperature in 
an intermediate region of 2R e < r < 10R Q (fig.l), hence, regulates the density scale height 
(or decreasing slope of density) in that region. The data based on LASCO/SOHO indicates 
that models considering larger r(=300s) are more likely. / max determines density decrease in 
a region of 1 ~ 2R Q . The data of CDS (< 1.2R Q ) and UVCS (~ 1.6i? Q ) exhibits the drastic 
decrease of density, which indicates that non-radial expansion is desired. Adjustment of the 
other parameters of the flow tube geometry (ri and a; eq.(14)) would give the still better 
fit. Although we do not further search the best parameter set to fit to the observation, the 
figure indicates that our model could reproduce the observed density profile by the choice of 
the appropriate parameters (-F Wi0 — (5 — 8) x 10 5 erg cm _2 s _1 , r m 300s, and / max 5). 



3.3.2. Temperature Distribution 

Figure 9 compares our results of temperature distribution with observation in the 
streamer region. Our models employ the same parameter sets as in fig. 8, and the data 
were electron temperature obtained from observations of a line ratio of Fe XIII/Fe X by 
the CDS and UVCS/SOHO (Parenti et al. 2000). Although the observed data are electron 
temperature, they are supposed to represent the plasma temperature because electron-ion 
equilibrium is attained in dense streamer region of r < 2R Q (Raymond et al. 1998). There- 
fore, it is reasonable to compare them to the results of our model considering one-fluid 
coronal plasma. A case employing (F Wi0 , r, / max ) = (7.8 x 10 5 ,300,5) gives reasonable peak 
temperature of ~ 1.5 x 10 6 K. However, none of our models can reproduce the observed loca- 
tion of T max . While the location is observationally inferred to be between 1.2i? Q and 1.6i?©, 
all of our models give r^ max < 1.2R & . This is because the dissipation length of the N- waves 
is essentially short, even though one considers long period- waves that are generated in the 
corona. We can summarize that acoustic waves excited in the corona could certainly heat the 
surrounding plasma to T > 10 6 K, however, they cannot maintain the high temperature till 
the sufficiently distant region by themselves. Therefore, the cooperations with other heating 
sources with larger dissipation length are necessary to explain the observed solar corona. 
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3.3.3. Velocity Distribution 

In fig. 10, we show the results of velocity distribution of the solar wind, with observational 
results in the low-latitude streamer (shaded). The observational data are from Sheeley et al. 
(1997), who determined velocity profile between 2 and 30 R & from measurements of about 65 
moving objects in the streamer belt. They used two different technic in deriving the results, 
whereas the shaded area displayed in fig. 10 are based on the straight-line fit method (an 
upper panel of fig. 6 in Sheeley et al. (1997); The shaded region is traced from that figure). 
Figure 10 indicates that the resultant velocity of our models is lower than the observed data 
in the whole region. The observation exhibits rapid acceleration in 3 ~ 5R & , while our 
results show gradual acceleration by thermal gas pressure, since the N-waves are damped in 
< 1.5i? Q and wave pressure cannot contribute to the acceleration of the wind flow, as shown 
in figs.l and 2. Other mechanisms are also required in acceleration of the low-speed wind. 

4. Summary and Discussions 

We have investigated the process of acoustic waves generated in the corona as a heating 
source especially in the mid- to low-latitude corona where the low-speed winds come from. 
We have found that the acoustic waves with r > 60s and F W) o > 2 x 10 5 erg cm _2 s _1 could 
certainly heat up the ambient plasma to T > 10 6 K by the dissipation of the N-waves even 
in a flow tube with the areal expansion of / max = 5, by balancing with the losses of the 
radiative cooling, the downward thermal conduction, and the adiabatic loss by the solar 
wind. Due to its dissipative character, the dissipation of N-waves effectively works in the 
inner corona and reproduces the density profile observed in the streamer region. However, it 
cannot contribute to the heating of the outer corona, since most of the wave energy is damped 
within a region of few-tenth of the solar radius even if one considers the waves with a long 
period of r = 300s. As a result, it is impossible to explain the observed temperature profile 
and flow velocity of the low-speed wind only by that process. Therefore, other mechanisms 
with the larger dissipation length should play a role cooperatively in the coronal heating and 
the acceleration of the low-speed wind. 

Small-scale reconnection events and spicules also predict generation of fast shock waves 
(Lee & Wu 2000; Lee 2001; Hollweg 1982), though we have concentrated on the role of the 
N-waves (slow shocks along the magnetic field line) excited by those events in this paper. 
The dissipation length of the fast shocks must be larger than N-waves with the identical 
period, since their phase speed is > va, which is much larger than that of the N-waves 
(~ c s ) for low-/? coronal plasma. Hence, cooperation with the fast shocks would give a more 
distant location of T max to match the observed temperature profile (fig.9), and would let the 
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acceleration of the wind flow continue in the outer corona, as suggested by observation of 
the streamer belt (fig. 10). 

Finally, we had better remark an issue on anisotropic and selective heating of ions. 
According to recent observations (Strachan et al. 2002), O VI ions in the streamer have high 
perpendicular kinetic temperature, though this is not so extreme as that observed in the 
high-speed wind. Our results show that the deposition of energy and momentum from the 
N-waves to the ambient gas is completed in the region of r < 1.5i? , where the particles are 
thermally well coupled. As a result, whatever anisotropies obtained by the N-wave heating 
would be wiped out and all the ions would have isotropic kinetic temperature. The process 
of the N-waves cannot explain the observed anisotropies, which also indicates that the other 
mechanisms causing the anisotropic heating have to operate simultaneously. 

We thank K. Shibata, K. Ohki, K. Omukai, K. Tomisaka, H. Saio, S. Nitta, S. Inutsuka, 
T. Kudoh, Y. Yoshii, T. Kajino, and A. Tohsaki as well as members of DTAP in NAOJ for 
many valuable and critical comments and Y. Mclean for improvement of presentation in this 
paper. The author is supported by the JSPS Research Fellowship for Young Scientists, grant 
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Fig. 1. — Variation of coronal wind structures on different sets of wave periods (in second) 
and the geometrical expansion factor, (r, / max ) = (300,5) (solid), (60,5) (dashed), and (300,1) 
(dotted) for the same F Wj0 = 7.8 x 10 5 erg cm _2 s -1 . On the left, distributions of temperature, 
flow velocity, and density from the inner boundary to 1AU (= 215i? ) are displayed from 
the top to the bottom. On the right, wave amplitude, a w , scaled by ambient sound velocity, 
Dissipation length, l w , of the waves, and heating per mass from the inner boundary to 3R & 
are presented from the top to the bottom. 
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Fig. 2. — Variation of coronal wind structures on three different F W; o = (3.2, 5.9, 10) x 10 5 erg 
cm _2 s _1 (solid, dashed, and dotted lines, respectively) for r = 120s and / max = 5. Each 
panel is same as fig. 1. 
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Fig. 3. — Relation between F Wi0 and T max on different sets of (r, f max ). 
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Fig. 5. — Relation between F w0 and (n p v )iau on /max = 1 and 5. Solid lines are results 
for cases adopting r = 300s, dotted lines, 120s, and dashed lines, 60s. Shaded region is the 
observational constraint by Withbroe (1988) for 'quiet corona'. 
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Fig. 6. — Fraction of three types of coronal energy loss, downward thermal conduction, 
radiative escape, and flow of the solar wind (mass loss), as a function of F w for the r = 120s 
waves in the flow tube with / max = 5. The non-labeled region between 'radiative escape' and 
'flow' indicates energy carried away to the outside of r out = 300-R© by outward conduction. 
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Fig. 7. — Exchange of energy among four components of energy flux, wave, radiation, conduc- 
tion, and flow (consult text for detail). The 'flow' component is increased by 10 times. Models 
adopt the same inputs of r = 210s and / max = 5, but different F w = (3.2, 5.9, 10) x 10 5 erg 
cm _2 s _1 (solid, dashed, and dotted lines, respectively) 
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Fig. 8. — Comparison of electron density distributions. Lines are model results and points 
are observational data. Models adopt the same F w = 7.8 x 10 5 erg cm _2 s _1 and different sets 
of (r(s),/ max ) = (300,5) (solid), (60,5) (dashed), and (300, f) (dotted). Pluses are data from 
CDS/SOHO observation of the mid-latitude streamer, asterisks from CDS/SOHO observa- 
tion of the equatorial streamer, diamonds from UVCS/SOHO observation of the mid-latitude 
streamer, and triangles from UVCS/SOHO observation of the equatorial streamer (Parenti 
et al. 2000). Squares are derived from observation of total brightness by LASCO/SOHO 
(Hayes, Vourlidas, & Howard 2001). 
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Fig. 9. — Comparison of electron density distributions. Respective lines represent results of 
the models adopting the same parameter sets as in fig. 8. Pluses are data from CDS/SOHO 
observation of the mid-latitude streamer, asterisks from CDS/SOHO observation of the 
equatorial streamer, diamonds from UVCS/SOHO observation of the mid-latitude streamer, 
and triangles from UVCS/SOHO observation of the equatorial streamer (Parenti et al. 2000). 
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Fig. 10. — Comparison of velocity profile of the solar wind. Respective lines represent results 
of the models adopting the same parameter sets as in fig. 8. Shaded region is observational 
data derived from measurements of about 65 moving objects in the streamer belt (Sheeley 
et al. 1997) (See text). 



